Structure of a prereaction complex between the nerve agent sarin, its biological target acetylcholinesterase, and the antidote HI-6.

Organophosphorus nerve agents interfere with cholinergic signaling by covalently binding to the active site of the enzyme acetylcholinesterase (AChE). This inhibition causes an accumulation of the neurotransmitter acetylcholine, potentially leading to overstimulation of the nervous system and death. Current treatments include the use of antidotes that promote the release of functional AChE by an unknown reactivation mechanism. We have used diffusion trap cryocrystallography and density functional theory (DFT) calculations to determine and analyze prereaction conformers of the nerve agent antidote HI-6 in complex with Mus musculus AChE covalently inhibited by the nerve agent sarin. These analyses reveal previously unknown conformations of the system and suggest that the cleavage of the covalent enzyme-sarin bond is preceded by a conformational change in the sarin adduct itself. Together with data from the reactivation kinetics, this alternate conformation suggests a key interaction between Glu202 and the O-isopropyl moiety of sarin. Moreover, solvent kinetic isotope effect experiments using deuterium oxide reveal that the reactivation mechanism features an isotope-sensitive step. These findings provide insights into the reactivation mechanism and provide a starting point for the development of improved antidotes. The work also illustrates how DFT calculations can guide the interpretation, analysis, and validation of crystallographic data for challenging reactive systems with complex conformational dynamics.

The Nature of Activated Non-classical Hydrogen Bonds: A Case Study on Acetylcholinesterase-Ligand Complexes.

Molecular recognition events in biological systems are driven by non-covalent interactions between interacting species. Here, we have studied hydrogen bonds of the CH⋅⋅⋅Y type involving electron-deficient CH donors using dispersion-corrected density functional theory (DFT) calculations applied to acetylcholinesterase-ligand complexes. The strengths of CH⋅⋅⋅Y interactions activated by a proximal cation were considerably strong; comparable to or greater than those of classical hydrogen bonds. Significant differences in the energetic components compared to classical hydrogen bonds and non-activated CH⋅⋅⋅Y interactions were observed. Comparison between DFT and molecular mechanics calculations showed that common force fields could not reproduce the interaction energy values of the studied hydrogen bonds. The presented results highlight the importance of considering CH⋅⋅⋅Y interactions when analysing protein-ligand complexes, call for a review of current force fields, and opens up possibilities for the development of improved design tools for drug discovery.

Hydroxyethylene isosteres introduced in type II collagen fragments substantially alter the structure and dynamics of class II MHC A(q)/glycopeptide complexes.

Class II major histocompatibility complex (MHC) proteins are involved in initiation of immune responses to foreign antigens via presentation of peptides to receptors of CD4(+) T-cells. An analogous presentation of self-peptides may lead to autoimmune diseases, such as rheumatoid arthritis (RA). The glycopeptide fragment CII259-273, derived from type II collagen, is presented by A(q) MHCII molecules in the mouse and has a key role in development of collagen induced arthritis (CIA), a validated model for RA. We have introduced hydroxyethylene amide bond isosteres at the Ala(261)-Gly(262) position of CII259-273. Biological evaluation showed that A(q) binding and T cell recognition were dramatically reduced for the modified glycopeptides, although static models predicted similar binding modes as the native type II collagen fragment. Molecular dynamics (MD) simulations demonstrated that introduction of the hydroxyethylene isosteres disturbed the entire hydrogen bond network between the glycopeptides and A(q). As a consequence the hydroxyethylene isosteric glycopeptides were prone to dissociation from A(q) and unfolding of the ß1-helix. Thus, the isostere induced adjustment of the hydrogen bond network altered the structure and dynamics of A(q)/glycopeptide complexes leading to the loss of A(q) affinity and subsequent T cell response.

Benefits of statistical molecular design, covariance analysis, and reference models in QSAR: a case study on acetylcholinesterase.

Scientific disciplines such as medicinal- and environmental chemistry, pharmacology, and toxicology deal with the questions related to the effects small organic compounds exhort on biological targets and the compounds' physicochemical properties responsible for these effects. A common strategy in this endeavor is to establish structure-activity relationships (SARs). The aim of this work was to illustrate benefits of performing a statistical molecular design (SMD) and proper statistical analysis of the molecules' properties before SAR and quantitative structure-activity relationship (QSAR) analysis. Our SMD followed by synthesis yielded a set of inhibitors of the enzyme acetylcholinesterase (AChE) that had very few inherent dependencies between the substructures in the molecules. If such dependencies exist, they cause severe errors in SAR interpretation and predictions by QSAR-models, and leave a set of molecules less suitable for future decision-making. In our study, SAR- and QSAR models could show which molecular sub-structures and physicochemical features that were advantageous for the AChE inhibition. Finally, the QSAR model was used for the prediction of the inhibition of AChE by an external prediction set of molecules. The accuracy of these predictions was asserted by statistical significance tests and by comparisons to simple but relevant reference models.

Divergent structure-activity relationships of structurally similar acetylcholinesterase inhibitors.

The molecular interactions between the enzyme acetylcholinesterase (AChE) and two compound classes consisting of N-[2-(diethylamino)ethyl]benzenesulfonamides and N-[2-(diethylamino)ethyl]benzenemethanesulfonamides have been investigated using organic synthesis, enzymatic assays, X-ray crystallography, and thermodynamic profiling. The inhibitors' aromatic properties were varied to establish structure-activity relationships (SAR) between the inhibitors and the peripheral anionic site (PAS) of AChE. The two structurally similar compound classes proved to have distinctly divergent SARs in terms of their inhibition capacity of AChE. Eight X-ray structures revealed that the two sets have different conformations in PAS. Furthermore, thermodynamic profiles of the binding between compounds and AChE revealed class-dependent differences of the entropy/enthalpy contributions to the free energy of binding. Further development of the entropy-favored compound class resulted in the synthesis of the most potent inhibitor and an extension beyond the established SARs. The divergent SARs will be utilized to develop reversible inhibitors of AChE into reactivators of nerve agent-inhibited AChE.

Targeting acetylcholinesterase: identification of chemical leads by high throughput screening, structure determination and molecular modeling.

Acetylcholinesterase (AChE) is an essential enzyme that terminates cholinergic transmission by rapid hydrolysis of the neurotransmitter acetylcholine. Compounds inhibiting this enzyme can be used (inter alia) to treat cholinergic deficiencies (e.g. in Alzheimer's disease), but may also act as dangerous toxins (e.g. nerve agents such as sarin). Treatment of nerve agent poisoning involves use of antidotes, small molecules capable of reactivating AChE. We have screened a collection of organic molecules to assess their ability to inhibit the enzymatic activity of AChE, aiming to find lead compounds for further optimization leading to drugs with increased efficacy and/or decreased side effects. 124 inhibitors were discovered, with considerable chemical diversity regarding size, polarity, flexibility and charge distribution. An extensive structure determination campaign resulted in a set of crystal structures of protein-ligand complexes. Overall, the ligands have substantial interactions with the peripheral anionic site of AChE, and the majority form additional interactions with the catalytic site (CAS). Reproduction of the bioactive conformation of six of the ligands using molecular docking simulations required modification of the default parameter settings of the docking software. The results show that docking-assisted structure-based design of AChE inhibitors is challenging and requires crystallographic support to obtain reliable results, at least with currently available software. The complex formed between C5685 and Mus musculus AChE (C5685•mAChE) is a representative structure for the general binding mode of the determined structures. The CAS binding part of C5685 could not be structurally determined due to a disordered electron density map and the developed docking protocol was used to predict the binding modes of this part of the molecule. We believe that chemical modifications of our discovered inhibitors, biochemical and biophysical characterization, crystallography and computational chemistry provide a route to novel AChE inhibitors and reactivators.